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ABSTRACT 

We present a catalog of ^14,000 X-ray sources observed by the ACIS instrument on the Chandra X-ray 
Observatory within a 1.42 square degree survey of the Great Nebula in Carina, known as the Chandra Carina 
Complex Project (CCCP). This study appears in a Special Issue of the ApJS devoted to the CCCP. Here, we 
describe the data reduction and analysis procedures performed on the X-ray observations, including calibration 
and cleaning of the X-ray event data, point source detection, and source extraction. The catalog appears to 
be complete across most of the field to an absorption-corrected total-band luminosity of ^lO'^'^'^ 
erg s^^ for a typical low- mass pre- main sequence star. Counterparts to the X-ray sources are identified in 
a variety of visual, near-infrared, and mid-infrared surveys. The X-ray and infrared source properties presented 
here form the basis of many CCCP studies of the young stellar populations in Carina. 

Subject headings: ISM: individual objects (Carina Nebula) — open clusters and associations: individual (Trl4, 
Trl5, Trl6) — stars: formation — stars: pre-main sequence — X-Rays: stars 



INTRODUCTION 



The Chandra Carina Complex Project (CCCP) is a 1.42 square degree survey of the Great Nebula in Carina 
with the Imaging array of the Advanced CCD Imaging Spectrometer (ACIS-I, Garmire et al.|[2003 ) on the Chandra 
X-ray Observatory. Townsley et al. (2011a) discuss the survey motivation and design, provide an observing log for 



the ACIS observations, and present several figures that place the footprint of the survey in an astronomical context. 
That work and the material presented here initiate a Special Issue of the ApJS devoted to the CCCP. 

We describe here the data reduction and analysis procedures underlying various data products that are widely 
used by the remaining studies in this series. Section [2] discusses the calibration and cleaning of the X-ray event 
data, point source detection, and source extraction. Properties of the resulting X-ray point sources are presented in 
Section [sj Visual and infrared (IR) counterparts to the X-ray sources are reported in Section |4j and the reliability 
of the counterpart identification is discussed in Section [sj Section [6] discusses the completeness and reliability of the 
X-ray catalog. 



DATA REDUCTION AND EXTRACTION METHODS 



The CCCP survey consists of 38 ACIS-I observations totaling 1.2 Ms organized into 22 pointings (Townsley 



et al. 2011a Table 1 and Figure 4). Chandra surveys of similar complexity include a ^2 Ms observation of the 
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Galactic Center region (Muno et al. 2009), the CHAMP survey of 130 extragalactic fields (Kim et al. 20071, the 
^1.8 Ms observation of the extragalactic C-COSMOS region ( Puccetti et al.|[2009 ), and a ^2 Ms observation of the 



extragalactic CDF-S field ( |Luo et al.|[2008l ) 



Our data analysis techniques are discussed at length by Broos et al. (2010), hereafter referred to as BIO; we 
briefly review the procedure here. Chandra-AGIS event data are calibrated and cleaned as described in Section 3 



of BIO Those procedures seek to improve the accuracy of event properties (individual event positions, alignment 
of the Chandra coordinate system to an astrometric reference, energy calibration) and seek to discard events that 
are likely caused by various instrumental background components. A cosmic ray artifact known as "afterglow']^ 
creates a particularly troublesome type of background — a group of events appearing at nearly the same location on 
the detector in nearly consecutive CCD frames. Since no single currently-available method for identifying afterglow 
events is appropriate for both weak and bright sources, we adopt a bifurcated workflow in which heavily-cleaned data 



are used for source detection and lightly-cleaned data are used for source extraction (see BIO Section 3 and Figure 
!)• 

Data from two detectors lying far off-axis in the ACIS S-array, available for most observations, were calibrated 
and cleaned as above. However, S-array data were not used for source detection and were not routinely extracted, 
due to the very poor angular resolution of the Chandra mirrors at large off-axis angles. In a few cases, very bright 



sources produced S-array data that were very useful for studies of individual stars (Townsley et al. 2011a Parkin et 



aI][20lT ) 



Candidate point sources are identified in each of the CCCP pointings individually using two methods (see B10[ 
Section 4.2). First, twelve images (energy bands 0.5-2 keV, 2-7 keV, 0.5-7 keV in combination with four pixel sizes) 
are searched by the standard Chandra detection tool, wavdetec^ ( [Freeman et al.|[2002 ) . Second, Lucy- Richardson 
image reconstruction (Lucy 1974) is performed on overlapping images (0.5 8 keV) that tile each 
CCCP pointing, sized (1.5' x 1.5') so that the Chandra PSF is relatively constant within each tile; 
peaks in the reconstructed tiles are adopted as candidate sources. Candidate sources are not obtained 
from observations in other bands (e.g., visual or IR). 

These candidate sources are extracted from all the observations in which they appear on the ACIS I-array using 
the ACIS Extracf\{AE) package (see |B10| Section 5). Sixty percent of CCCP sources were observed multiple (2-7) 
times, as illustrated by the CCCP exposure map, which is shown in Figure [8] here. Figure 4 in [Townsley et al.| 



(2011a), and Figure 2 in BIO The range of off-axis angles, and thus the range of extraction aperture sizes, found 
among the multiple observations of a single source is often small — when observations share the same pointing or 
when observations barely overlap — but can be large when observations moderately overlap. 

The positions of source candidates are updated with AE estimates, using a subset of each source's extractions 



chosen to minimize the position uncertainty (see BIO Section 6.2 and 7.1). A source significance statistic is calculated 
with respect to the local background level (see B10[ Section 4.3), using a subset of each source's extractions chosen to 



maximize that significance (see BIO Section 6.2). This decision to ignore some observations of a source effectively 
prevents the very poor point spread function (PSF) of an off-axis observation from spoiling the detection or position 
estimate of a source that was also observed with a very good PSF on-axis. As the Carina field has bright astrophysical 
diffuse emission with complex structure, the use of a local rather than global background source existence criterion 
is important. Only candidates found to be significant are accepted as X-ray sources; expressed in terms of source 
properties reported in Table [ij a detection requires both SrcCounts.t > 3 and ProbNoSrc.min < 0.01. 
Iterative pruning of source candidates and re-extraction continues until no candidates are found to be insignificant 



(see BIO Section 4.1). Our selection of the ProbNoSrc_min threshold is discussed further in Section 6.2 



Photometry, spectra, light curves, and a variety of apparent source properties (Section|3| are derived for sources 
in the final catalog, using a subset of each source's extractions chosen to balance the conflicting goals of minimizing 



photometric uncertainty and of avoiding photometric bias (see BIO Section 6.2) 



The CCCP survey region includes several rich young stellar clusters, and some sparse but compact stellar groups. 



^ http : 1 1 cue . h arvard . edu/clao/why/ afterglow ■ html 



http : //asc ■ harvard . edu/ciao/dognload/doc/ detect.manual 



http: //www. astro .psu. edu/xray/acis/acis_analysis .html 



where source crowding can be important. In such cases, the source detection is largely based on the maximum 
likelihood image reconstruction where the blurring effects of the PSF are reduced. Photons are then extracted in 
apertures based on the local PSF shape but reduced in size to avoid overlap with adjacent source extraction regions. 
While the procedures are well-established in |B10| with little subjective judgment involved, it is important to recognize 
that the source identifications and extractions in very crowded regions represent only one plausible interpretation of 
the data. 



The Chandra data analysis system, CI AO, (Fruscione et al. 2006), the SAOImage ds9 visualization tool (Joye 
fc: Mandel||2003 ), and the Interactive Data Languag^ (IDL) are used throughout our data analysis workflow, from 
data preparation through science analysis. The CCCP data were reduced with CIAO version 4.0.2, using CALDB 
version 3.4.2 for event processing and CALDB version 4.1.1 for construction of calibration data products. Models of 
the Chandra-AClS PSF were constructed (see BIO Appendix C) using version 4.3 of the MARX mirror simulatoiF] 



3. APPARENT X-RAY POINT SOURCE PROPERTIES 



A total of 14,36cj3 point sources were identified and extracted. The CCCP observations are not appropriate for 
studying the luminous blue variable 77 Carinae, and it is omitted from our catalog even though it is a strong X-ray 
source; Townsley et al. (2011a) provide references to more suitable observations of this remarkable star. The CCCP 



catalog is depicted in Townsley et al. (2011a Figure 3 and 4), and in Broos et al. (2011 Figure 4). No single figure 



can adequately represent such a wide (1.4 square degrees) Chandra field since sources as close as <1" separation can 
be detected. However, the electronic versions of those figures can reveal, when zoomed, detail not visible in most 
printed versions. 



Estimates for many apparent (not corrected for absorption) properties of the X-ray sources (see BIO Section 7) 



are provided in a table published electronically and available at Vizier ( Ochsenbein et al.]|2000 ). Since the table has 
many columns, a stub cannot be conveniently shown in print. Instead, column names that are appropriate for Vizier 
and column descriptions are listed in Table [l] The suffixes "_t", "_s", and "_h" on names of photometric quantities 
designate the total (0.5-8 keV), soft (0.5-2 keV), and hard (2-8 keV) energy bands. The SrcCounts and NctComits 
quantities characterize the extraction; correction for finite extraction apertures is applied to the ancillary reference 



file (ARF) calibration products (see BIO Section 5.3). Thus, the lowest-level calibrated photometric quantity that 



can be used to compare sources would be apparent photon fiux (see |B10[ Section 7.4) 

^photon = NetCounts/McanEffcctivcArea/ExposureTimeNominal 
which has units of photon cm^^ s~^. Rough, model-independent estimates for apparent energy flux are provided as 



(1) 



EnergyFlux = 1.602 x 10 ^ MedianEnergy x Fphoto 



(2) 



in units of erg cm s , where the constant 1.602 x 10 arises from the conversion between keV and erg ( Getman et 



al.|2010 ). Table notes provide additional information regarding the definition of source properties. Model-dependent 
estimates for absorption and intrinsic energy flux are provided by Broos et al. (2011 Table 8) for sources that are 
likely to be low-mass young stars in the Carina complex. 



http: //www. ittvls . com/idl 



http : //space . mit . edu/cxc/marx/ 



The catalog used in other CCCP studies contains 14,368 sources; one additional source (the last reported in Table [ijl was identified 
too late in the data analysis process to be included in those studies. 



Table 1. CCCP X-ray Sources and Properties 



Column Label Units Description 



Name 

Label^ 
RAdcg 
DEdcg 
PosErr 
PosTypc 



deg 
deg 
arcsec 



lAU source name; prefix is CXOGNC J 

{Chandra X-ray Observatory Great Nebula in Carina) 

source name used within the CCCP project 

right ascension (J2000) 

dechnation (J2000) 

1-a error circle around ( RAdcg, DEd eg) 

algorithm used to estimate position JBIO Section 7.1) 



ProbNoSrc_min 
ProbNoSrc_t 
ProbNoSrcs 
ProbNoSrc_h 



smallest of ProbNoSrc_t, ProbNoSrcs, ProbNoSrc_h 



p-value for no-source hypothesis i BIO Section 4.3) 
p-value for no-source hypothesis 
p-valuc for no-source hypothesis 



ProbKS-single'^ 



ProbKS_merge'^ 



smallest p-value for the one-sample Kolmogorov-Smirnov 
statistic under the no- variability null hypothesis within a 
single-observation 

smallest p-value for the one-sample Kolmogorov-Smirnov 
statistic under the no-variability null hypothesis over merged 
observations 



ExposuroTimeNominal 
ExposureFr action'^ 
NumObservations 
NumMcrged 

McrgcBias 



total exposure time in merged observations 

fraction of ExposurcTimeNominal that source was observed 

total number of observations extracted 

number of observations merged to estimate photometry prop- 
erties 

fraction of exposure discarded in merge 



Theta.Lo 

Theta 

Thcta.Hi 



arcmm 
arcmin 
arcmin 



smallest off-axis angle for merged observations 
average off-axis angle for merged observations 
largest off-axis angle for merged observations 



PsfFraction ■ ■ ■ average PSF fraction (at 1.5 keV) for merged observations 

SrcArca (0.492 arcsec)^ average aperture area for merged observations 

AfterglowFraction^ ■ ■ - suspected afterglow fraction 



SrcCounts_t 
SrcCounts_s 
SrcCounts_h 



count 
count 
count 



observed counts in merged apertures 
observed counts in merged apertures 
observed counts in merged apertures 



BkgScaling 



scaling of the background extraction (IBlOl Section 5.4) 



BkgCounts_t 
BkgCounts_s 
BkgCounts.h 



count 
count 
count 



observed counts in merged background regions 
observed counts in merged background regions 
observed counts in merged background regions 



NetCounts.t 
NetCounts_s 
NetCounts.h 



count 
count 
count 



net counts in merged apertures 
net counts in merged apertures 
net counts in merged apertures 



MeanEffectivcArea.t* 
MeanEffectiveArca_s 
McanEffectivcArea_h 



cm" count photon" 
cm count photon" 
cm^ count photon" 



mean ARE value 
mean ARE value 
mean ARE value 



McdianEncrgy.t^ 

McdianEnergy_s 

McdianEnergy_h 



keV 
keV 
keV 



median energy, observed spectrum 
median energy, observed spectrum 
median energy, observed spectrum 



EnergyElux.t 
EnergyElux_s 
Energy Elux_h 

NetCounts_Lo_t^^ 
NetCounts_Hi_t 



erg cm 
erg cm" 
erg cm" 

count 
count 



max( Ene rgy Elux_s,0) -|- max(EnergyElux_h,0) 
See Eqn. 2 
See Eqn.p 



1-sigma lower bound on NetCounts.t 
1-sigma upper bound on NetCounts_t 



NetCounts_Lo_s 
NetCounts_Hi_s 



count 
count 



1-sigma lower bound on NetCounts_s 
1-sigma upper bound on NetCounts_s 



NetCounts_Lo_h 
NetCounts_Hi_h 



count 
count 



1-sigma lower bound on NetCounts_h 
1-sigma upper bound on NetCounts_h 



Note. — Column labels are shown as they appear in Vizier, where formatted labels are not possible. 

Note. — The first 14,368 sources — the catalog used in other CCCP studies — are sorted by RA; the 14,369*'* source 
was identified too late in the data analysis process to be conveniently inserted at its rightful position in the table. 



Note. — The suffixes "_t", "_s" , and "_h" on names of photometric quantities designate the total (0.5—8 keV), soft 



(0.5—2 kcV), and hard (2—8 kcV) energy bands. 
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Note. — Source significance quantities (ProbNoSrc_t, ProbNoSrc_s, ProbNoSrc_h, ProbNoSrc_min) are computed 
using a subset of each source's extractions chosen to maximize significance ) |B10| Section 6.2). Source position 
quantities (RAdeg, DEdeg, PosErr) are computed using a subset of each source's extractions chosen to minimize the 
position uncertainty | |B10| Section 6.2 and 7.1). All other quantities are computed using a subset of each source's 
extractions chosen to balance the conflicting goals of minimizing photometric uncertainty and of avoiding photometric 
bias jB10| Section 6.2 and 7). 

^Source labels identify a CCCP pointing | |Townsley et ar||2011a] Table 1); they do not convey membership in 
astrophysical clusters. 

'^In statistical hypothesis testing, the p-value is the probability of obtaining a test statistic at least as extreme as 
the one that was actually observed when the null hypothesis is true. 

'^See |Bl'o| Section 7.6 for a description of the variability metrics, and caveats regarding possible spurious indications 
of variability using the ProbKS_merge metric. 

'^Due to dithering over inactive portions of the focal plane, a Chandra source is often not observed during some 
fraction of the nominal exposure time. (See http://cxc.harvard.edu/ciao/why/dither.htmll) The reported quantity 
is FRACEXPO produced by the CI AO tool mkarf. 

^ Since the extracted event data are lightly-cleaned to avoid removing legitimate X-ray events from bright sources 
||B10l Section 3), some background events arising from an effect known as "afterglow" ( http : //cxc . harvard, edu/ciao/l 
why/af terglow . html } will remain. We attempt to identify afterglow events using the tool ae_af terglow_report (see 
the AE manual), and report the fraction of extracted events attributed to afterglow. 

^The ancillary response file (ARF) in ACIS data analysis represent both the effective area of the observatory and 
the fraction of the observation for which data were actually collected for the source (ExposureFraction). 

^MedianEnergy is the quantity ENERG_PCT50_OBSERVED, the median energy of extracted events, corrected 
for background l |B10[ Section 7.3). 

'^Confidence intervals (68%) for NetCounts quantities are estimated by the CI AO tool aprates l |http : //asc. harvard. | 
edu/ciao/ahelp/aprates .html ) . 
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High-quality spectra are available for relatively few CCCP sources. Table [2] shows the number of CCCP sources 
that have spectra with various levels of quality (characterized by ranges of net counts in the total band, NetCounts_t) ; 
the tallies are reported separately for a set of very significant "primary" sources, defined as having a low probability 
of being a spurious source (ProbNoSrc_min in Table [l] less than 0.003) and the complementary set of less significant 
"tentative" sources (0.003 < ProbNoSrcmin < 0.01). 

The 40 sources with the highest-quality spectra (more than 500 net counts) are listed in Table [s] Most are 
known or candidate OB stars that are discussed in other CCCP studies ( Naze et al.|2011 Povich et al.|2011a[ Parkin 



et al.|2011[ ). The three known Wolf-Rayet stars in Carina (WR 25, WR 24, WR 22) are discussed by |Townsley et al 



( |2011a[ ). The neutron star 104608.71-594306.4 (El_85) is discussed by |Hamaguchi et al.| pOOO] ). The 12 remaining 



sources in the table exhibit the time variability and high-temperature spectra (shown in Figure [T]) characteristic of 
pre-main sequence stars during flares attributed to magnetic connection events; their inferred time-averaged intrinsic 
luminosities range from ^lO'^'^'^ erg s^^ to ~'10'^^ '' erg s^^ in the total-band (0.5-8 keV). 

Photon-counting detectors, such as ACIS, can suffer from a non-linearity known as photon pile-u;^wh.en multiple 
X-ray photons arrive with a separation in time and space that is too small to allow each to be detected as a separate 
X-ray event. Pile-up effects include photometric and spectroscopic miscalibration of the observation. Total-band 
photometry is underestimated, because multiple photons interact to produce only one or zero events. The shape 
of the detected X-ray spectrum is hardened, because the energy assigned to a piled event will represent that from 
multiple photons. 

AE screens all point source extractions for photon pile-up by estimating the observed count rate falling on an 
event detection cell of size 3x3 ACIS pixels, centered on the source position. Pile-up was confirmed and quantified for 
the at-risk extractions so identified, using an experimental Monte Carlo forward-modeling approach that reconstructs 
a pile- up free ACIS spectrum from a piled ACIS observation. Appendix |^ describes this method and illustrates its 
application to CCCP sources. Table [4] lists extractions of CCCP sources — all known massive stars — that were con- 
firmed to have significant pile-up using this method]^ Column (6) reports the ratio of the pile-up free to observed 
(piled) count rate in the total (0.5-8 keV) energy bandj^ For these piled sources. Table [l] reports total-band pho- 
tometry (columns SrcCounts_t and NetCounts_t) as observed, (without correction) and all other quantities that are 
affected by pile-up are omitted. 



http : //cxc . harvard . edu/ ciao/why/pileup_intro . html 

Although source 104544. 13-592428. 1=C4_2144=HD 93403 had a very high event rate (0.17 count s l), reconstruc- 
tion of the spectrum indicated very little pile-up (a value of 1.02 in column 7 of Table [l]| because the source was 
observed very far off-axis (8') and thus had only a modest event rate per detection cell. 

® We choose not to use the terms "pile-up fraction" or "pile-up percentage" because the ACIS community has several conflicting 
definitions for those terms; see Section 1.2 in The Chandra ABC Guide to Pileup ) |http: //cxc. harvard ■ edu/clao/do«nload/doc/pileup_| 
abc.pdf ). 



Table 2. Quality of Spectra for CCCP Sources 



NetCounts.t Number of sources 
Primary Tentative 

(count) 

(1) (2) (3) 



<10 


6852 


1424 


10- 50 


4934 


70 


50- 100 


677 





100- 500 


371 





500-1000 


24 





>1000 


16 
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Table 3. CCCP sources with >500 net counts in I-array extractions 



CXOGNC .1^ 


Labol"- 


NetCounta.t 
(count) 


Name 


SpTypc 


CCCP citation 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 



104410 


39 


694311 


1 


CTrl6-280 


24609 


WR 25 (HD 93162) 


WN6h + OB? 


104544 


13 


692428 


1 


C4_2144 


19161 


HD 93403 


05.5(t) 


104445 


04 


693354 


6 


CTrl4_3535 


9274 


HD 93250 


03.5V((f+)) 


104422 


91 


695935 


9 


C2_llll 


4169 


QZ Car (HD 93206) 


09.5I 


104357 


47 


593251 


3 


CTrl4_1925 


2977 


HD 93129A 


02It» 


104352 


25 


600704 





C2_547 


2664 


WR 24 (HD 93131) 


WN6-A 


104508 


23 


694607 





CTrl6-3102 


1909 


CI* Trumplcr 16 MJ 496 


08.5V 


103909 


94 


694714 


5 


SB2_6 


1864 




candidate OB star 


104505 


90 


594006 





CTrl6-2805 


1663 


HD 303308 


04V((t+)) 


104752 


54 


600215 


2 


C6_1260 


1673 






104456 


27 


693830 


4 


CTrl6_1450 


1468 






104433 


74 


694415 


4 


CTrl6-749 


1408 


HD 93205 


03.5V((t+)) + 08V 


104246 


63 


601207 





SP2_103 


1403 




candidate OB star 


104402 


75 


693946 





CTrl6-172 


1278 




candidate OB star 


104154 


91 


694123 


6 


Cl_99 


1249 




candidate OB star 


104815 


18 


694319 


7 


E4_67 


1084 






104243 


76 


693954 


2 


Cl-549 


976 






104441 


SO 


694656 


4 


CTrl6_1028 


976 


CI Trumplcr 16 100 


05.5V 


104055 


32 


694239 


7 


SB1_81 


965 






104417 


64 


695350 


6 


C2_995 


822 






104436 


23 


600529 





C2_1238 


820 


HD 93222 


08III((t)) 


104356 


16 


693250 


9 


CTrl4_1731 


716 






104357 


65 


693253 


7 


CTrl4_1963 


685 


HD 93129B 


03.5V((f+)) 


104609 


90 


602635 


5 


SP4_114 


680 






104837 


74 


601325 


7 


SP3_222 


668 


LS 1932 (HD 93843) 


05V 


104117 


60 


694037 





SB1_176 


664 


WR 22 (HD 92740) 


WN7-A 


104516 


62 


694337 


1 


CTrl6-3341 


624 


CI Trumplcr 16 112 


05.5/6V((f+?p)) + B2III/V 


104512 


86 


694942 


3 


CTrl6_3272 


619 






104455 


79 


695826 


9 


C3-210 


618 






104343 


29 


602309 


5 


SP2_52S 


600 






104457 


51 


695429 


5 


C3_225 


674 




candidate OB star 


104605 


70 


695049 


5 


El_74 


644 


LS 1886 (HD 305525) 


04V 


104615 


77 


692728 





C6_759 


629 






104408 


84 


693434 


4 


CTrl4_2688 


622 


HD 93161a 


08V + 09V 


104608 


71 


694306 


4 


Bl_85 


522 






104354 


40 


693257 


4 


CTrl4_1506 


519 


HD 93128 


03.5V((f+)) 


104220 


83 


690908 


6 


C7_58 


511 




candidate OB star 


104712 


63 


600550 


8 


C6_846 


502 


HD 93632 


05111(f) 


104557 


13 


595643 


1 


C3_1477 


600 


Hen 3- 485 





t al. 20lla1) 
ct al- 2011) 
ct al, ■2ULL) 
ct al, 2()ll;| 
ct al- ■JULLi 
t al. 2()llaJ 
ct al- 20111) 
t al. 2Ull£ii 




^Source labels identify a CCCP pointing jTownslcy et al.|2011a| Tabic 1); they . 
^lAV source name; prefix is CXOGNC J {Chandra X-ray Observatory Great Nebula 



CarinaJ 



L astrophysical clusters. 



Table 4. Sources exhibiting photon pile-up 



CXOGNC J 


Label 


Name 


Observation 


e 


PsfFraction 


Correction 


(1) 


(2) 


(3) 


(4) 


{') 
(5) 


(6) 


(7) 


104357.47-593251.3 


CTrl4_1925 


HD 93129A 


4495 


0.3 


0.61 


1.27 


104410.39-594311.1 


CTrl6_280 


WR 25 (HD 93162) 


6402 


4.8 


0.90 


1.25 


104433.74-594415.4 


CTrl6_749 


HD 93205 


6402 


1.8 


0.50 


1.16 


104422.91-595935.9 


C2_llll 


QZ Car (HD 93206) 


9482 


3.4 


0.90 


1.11 


104445.04-593354.6 


CTrl4_3535 


HD 93250 


4495 


6.3 


0.89 


1.05 



— Col. (1): lAU designation {Name in TablejT} 
Source name used within the CCCP project (Label in Table Q 
Object common name 
Observation ID 

Off-axis ang le {Theta in Table [T]) 

Fraction of the PSF (at 1.497 keV) enclosed within the extraction region (PsfFraction in Table A 
reduced PSF fraction (significantly below 90%) indicates that the source is in a crowded region. 
Col. (7): Estimated ratio of pile-up free to observed (piled) count rates in the 0.5—8 keV energy band. 
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Fig. 1. — High-quality spectra for sources from Table [3] expected to be pre-main sequence stars. The spectra are 
modeled as one-tempcrature or two-temperature thermal plasmas with abundances frozen at the values adopted by 
the XEST study ( Giidel et al.||2007[ ), relative to [Anders fc Grevessej ( |1989[ ), scaled to [Wilms et al.| ( |2000[ ), using the 
tbabs absorption code. The model is implemented in XSPEC as tbabs*(apec) or tbabs* (apec+apec) . Column 
density and plasma temperatures are shown on each panel in units of 10^^ cm^^ and keV. 



COUNTERPARTS TO X-RAY POINT SOURCES 



Many studies in this Special Issue rely on counterparts to CCCP sources observed in other bands. [Broos et al.[ 



(2011 ) assign to each source a probability of Carina membership based in part on near- and mid-infrared photometry 
and on visual spectroscopy. CCCP studies of massive stars rely on visual spectroscopy for conclusive identification 
Naze et aL][2011[ Parkin et al.[[20ir i and on visual or infrared photometry to identify candidates 



Evans et al. 2011 



Povich et al.[[20lTa| . Both X-ray and infrared data are essential to studies of pre-main sequence 



stars ( [Preibisch et~al][2011[ [Wang et al.[[2011[ [Wolk et al.[[2011 1 . 



Table[5|lists 19 counterpart catalogs that we have correlated with the CCCP source list. The Skiff, KR, PPMXL, 
UCAC3 (and its associated BSS), 2MASS, and Spitzer catalogs cover the entire CCCP field. The HAWK-I, SOFI, 
NACO, and Sana near-infrared catalogs cover portions of the field including the richest clusters. The CMD, DETWC, 
MDW, MJ, CP, and DAY visual photometry catalogs survey small fields on the major clusters. The X-ray study of 



Tr 16 by Albacete-Colombo et al. (2008) is discussed by Wolk et al. (2011 1 



The CCCP catalog was matched to counterpart catalogs one at a time using a simple algorithm, described in 
BIO (Section 8), in which the maximum acceptable separation between an X-ray source and a counterpart is based 
on the individual source position errors assuming Gaussian distributions, scaled so that ~99% of true associations 
should be identified as matches. When multiple sources in the counterpart catalog satisfy the match criterion, the 
closest one is adopted. The performance of the algorithm is studied in Section [5] 

For each counterpart catalog in TablejSj Neat gives the approximate number of catalog entries lying in the CCCP 
field, and Ncccp gives the number of catalog entries that are successfully matched with X-ray sources. Preliminary 
matching runs were used to estimate and remove shifts (but not rotations) between each counterpart catalog and an 
appropriate astrometric reference, usually the Naval Observatory Merged Astrometric Datase1p°] (NOMAD, |Zacharias 



et al. 2004); the ACIS data have also been aligned to NOMAD. These catalog offsets are given in columns (6) and 
(7). Median separations among the identified counterparts are reported as rmodian in Table [sj For most of the 
catalogs, >80% of the counterparts have separations less than 1". Figure [2] maps the CCCP sources with identified 
counterparts for some of the catalogs listed in Table [5] 

Selected information about the identified counterparts in all the catalogs in Table [5] are reported in a single 
table with 14,368 rows representing the CCCP sources. The table is available in the electronic version of this paper, 
and at Vizier ( Ochsenbein et al.|2000 ). Column names and descriptions are listed in Table [6j The first four columns 
are reproduced from Table [T] to identify each CCCP source. The Identifier, SpType, and SpRef columns provide an 
historical name for the source, report the spectral type we have adopted in the CCCP, and report a reference for that 



spectral type. These quantities are taken from Gagne et al. (20111 when available, and from Skiff (2009) otherwise 



Columns named "Ofi'set(*)" report distances between CCCP and counterpart sources, after astrometric align- 
ment of the catalogs. The remaining columns appear in groups corresponding to each counterpart catalog. Within 
each group, a column named "Id(*)" reports a unique identifier in the counterpart catalog. For the convenience of 
the reader|^ additional columns give photometry or other information from the counterpart catalog}^ 



http: //www. nofs .navy .mil/nomad/ 



Ideally, Table [g] would report only the counterpart identifiers, and the reader would retrieve whatever set of previously published 
counterpart columns are desired at the time Table |6] is downloaded. However, no astronomy catalog server known to the authors can 
perform this "join" service. 

No screening of the re-published columns has been performed, e.g., removal of photometry estimates that are flagged as low quality 
by the original authors. Readers interested in counterpart properties for detailed studies are advised to identify CCCP sources in the 
original catalogs (using the "Id" entries in Table |6]l and interpret that counterpart data appropriately. 



Table 5. Summary of catalogs correlated with Chandra Carina sources 



Catalog 


Scope 


Reference 




Ncccp 


AR.A. 


ADcc 














(") 


(") 


(") 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 



Skiff Visual spectral types 

KR Visual photometry 

PPMXL CCD proper motions (PMs) 

UCAC3 CCD PMs 

BSS Bright star PMs 

CMD Photographic PMs, Tr 14, Tr 16, Cr 232 

DETWC Visual photometry, Tr 14 & 16 

MDW Visual spectral types, Cr 228 

MJ Visual photometry, Tr 14 & 16 

CP High-mass photometry, Cr 228 

DAY Low-mass photometry, Cr 228 

HAWK-I Deep near-infrared photometry 

2MASS Shallow near-infrared photometry 

SOFI Deep near-infrared photometry, Tr 14 

NACO Deep near-infrared photometry, Tr 14 

Sana Deep near-infrared photometry, Tr 14 

SpVela Mid-infrared photometry (Spitzer) 

SpSmith Mid-infrared photometry (Spitzer) 

AC ACIS observation of Tr 16 



ISkifFV2009l 



Kharcheii ko Roeserj j2009| 

Uoeser et aFl 20101 
Zacharias et al. (201ojl 
Urban et al. (2004) 
Cudworth et ai. ( 19 93| 



DcGioia-Eastwood e t~al.| | |200l| 

Massey et al. (2001)^ 

iMassey & Johnson (199 311 
Carraro & Patat ( 20011 



Dclgado ct 


al. 




HAWKl-1 


206S 




Cutri ct al 


(20 


m 


Asccnso ct 


al. 1 


2007 


Asccnso ct 


al. 1 


2007 



Sana ct al. (20101 
Povich ct al. (2011b) 

STmt h ct al. (20101 

Albacctc-Colombo et al. 



I 2008 I 



271 
363 
30286 
23557 
47 
577 
853 

73=^ 
768 
1112 
152 
595148 
130164 
4739 
178 
1955 
130871 
97901 
1035 



130 
114 
1006 
1471 
18 
141 
237 
35 
178 
27 
17 
6583 
6194 
849 
54 
371 
6543 
3811 
935 



0.00 
0.05 
0.00 
0.00 
0.00 
0.62 
0.53 

0.43 
-0.64 
-0.05 
-0.01 
0.00 
0.02 
-0.01 
-0.06 
0.00 
0.00 
0.19 



0.01 
0.02 
0.00 
0.00 
0.00 
0.59 
0.16 

0.27 
0.11 

-0.03 

-0.02 
0.00 
0.00 

-0.01 
0.04 
0.00 
0.00 

-0.01 



0.19 
0.22 
0.41 
0.33 
0.16 
0.33 
0.27 
0.19 
0.41 
0.68 
0.81 
0.31 
0.37 
0.18 
0.17 
0.16 
0.38 
0.44 
0.18 



'^Thc reported positions for two of the 75 MDW sources do not appear to coincide with any 2MASS source. We ignored those sources, and 
adopted 2MASS positions for the 73 remaining. 



- 10 - 



Table 6. Counterpart Properties 



Column Label 


Units 


Description 


1x1 !^ m 






T,ahpl 




Jt— <3r^nTT*fi nnmp ncfin As/itrnn trip l,i;l,T-' ■r\Tr»iPr*t 

j\. L Chy OKJUHjC IICLIIIC LIOCU WJUlllll lillC V_y ± ^JLUJKjKjIi 


WiiaCLl OKlli 1 


3.rCS6C 


offset from. m.atching ACIS source 


WlibCLl xvXv, ) 


arcsGC 


offset from matching ACIS source 




arcsGc 


offset from matching ACIS source 




arcsec 


offset from matching ACIS source 




arcs6C 


offset from matching ACIS source 




arcs6C 


offset from matching ACIS source 


o ff<?p t ( Ti F,T wn 1 


dl UbtJU 


UliotJt llUlIi lliatUlilllg .rVOiO bULilUtJ 


Wlibt;i;l^iVli_-' VV ) 


arcsGC 


offset from. m.atching ACIS source 




arcsec 


offset from matching ACIS source 




arcsec 


offset from matching ACIS source 




arcsec 


offset from matching ACIS source 






i^frGPt TTOvn TTi cit r"}! in cr AflT^ QrinTr*P 




arcsec 


offset from matching ACIS source 


Offset I'fiOFP 


Oil- ^o^^ 


nffspt froTTi Tnatfhi'nP' ACTS snnrpp 


Offset (NACO) 


arcsec 


offset from matching ACIS source 


Offset(Sana) 


arcsec 


offset from matching ACIS source 


Offset (SpVela) 


arcsec 


offset from matching ACIS source 


Offset (SpSmith) 


arcsec 


offset from matching ACIS source 


Offset(AC) 


arcsec 


offset from matching ACIS source 


Identifier 




historical star name 


SpType 




spectral type 


SpRcf 




reference for spectral type 


Id(KR) 




row identifier {j'ecno in published catalog) 


Bniag(KR) 


mag 


B band magnitude 


Vmag(KR) 


mag 


V band magnitude 


Id(PPMXL) 




row identifier {ipix in published catalog) 


Id(UCAC3) 




row identifier {^3UC in published catalog) 


Id{BSS) 




row identifier i^ID in published catalog) 


Id(CMD) 




row identifier [recTio in published catalog) 


Vinas(CMD) 




\^ Ijfiiid iiifi^iiitndc 


B-\'((.'MiJ) 




J_)-\ color 


IdrDRTWC"! 




row identifier ( RA I2000-f-DE 12000 in nublished cataloff^ 


Vmag(DETWC) 




V band magnitude 


U-B(DETWC) 


mag 


U-B color 


B-V(DETWC) 


mag 


B-V color 


Id(MDW) 




row identifier ( Ncityib in published catalog) 


TH I'M 11 




1 UW lUtJllbiXlCX \iJVUf 111 pLiUlibllCU UatcUlUg J 




mag 


\/ nanri m a crni "1" n ri P 

V UCIiXXVJ. IXXdi^XXXli uuc? 


U-B(MJ) 




IJ-B color 


B-V(MJ) 


map" 


B-V color 


Id(CP) 




vi^iAs? iH OTi("iTi rif 1 11) in T^nnliGriPri r*Jitali^cri 

i W \y 1*_XL>11L111^X IJ J_y IXX LI 1_1 LIXXoIXC^U K^GjvCKIKJp^ t 


Tlmap-COPI 


mag 


T T Q n vn QrmTi"ii/Hti 

\j uaiici iixcigixi t U-Cic 


Rmao-rOPI 

J_J lllClg I V—' A 1 


mag 


R nanH mn crnitnrip 

J— > I^CLXXU XXXdiglllU LiUC 




mag 


V band magnitude 




mag 


R band magnitude 




mag 


I band magnitude 


Id (DAY) 




TOW iHpTitifipr ( f^fiTfip in nnbliciViprl fatalop'l 

X \J w xu^xiuxxx^x yiyijbiiv^ xix vj L^xioxxcvx \^<xh<\iL\J^ I 


Vmag(DAY) 


mag 


Johnson V magnitude 


U-B(DAY) 


mag 


Johnson U-B colour index 


B-V(DAY) 


mag 


Johnson B-V colour index 


V-Rc(DAY) 


mag 


Johnson-Cousins V-Rc colour index 


V-Ic(DAY) 


mag 


Johnson-Cousins V-Ic colour index 


AV(DAY) 




Absorption in V band 


Id(HAWK-I) 




row identifier {RA-j-DEC in published catalog) 


Id(2MASS) 




row identifier {Designation in published catalog) 




Fig. 2. — CCCP matches to selected catalogs in Table [5] The electronic version of the figure can reveal, when 
zoomed, detail not visible in most printed versions. Coordinates here and for all subsequent images are celestial 
J2000. 
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4.1. Chandra detections of Carina stars with visual spectroscopy and photometry 



The General Catalog of Stellar Spectral Classifications ( "Skiff" in Table [5]) is an all-sky collection of stars with 
MK spectral types based on visual spectroscopy gathered from heterogeneous sources in the literature; 271 Skiff 



sources are found in the CCCP field. Gagne et al. (2011) and Naze et al. (2011) study 200 O and B stars selected 



from the Skiff list; for a few of these stars we adopted positions and/or spectral types different from those reported 
by |Skiff| ( p009| ( [Gagne et al.|2011| . Three Wolf-Rayet stars in the field are discussed by |Townsley et al.| ( |2011a] ) 



stars with spectral types A, F, or K are assumed to lie in the foreground ( Broos et al.||2011 ). A considerable fraction 
of the known OB stars in the region are not detected in the CCCP (Gagne et al. 2011) and a considerable number 



of likely OB members are discovered in the CCCP and do not yet have spectroscopic classifications (Povich et al. 



2011a). Thus, the Skiff-CCCP counterparts can not be viewed as well-defined or complete in any criteria. 



Kharchenko & Roeser ( 2009| ) (KR) provide an all-sky catalog of visual-band measurements compiled from the 
literature for more than 2.5 million stars brighter than V ~ 14 mag; we identify 114 of these in the CCCP catalog. 



Massey fc Johnson] ( |1993[ ) (MJ) and [DeCioia-Eastwood et al] ( |2001[ ) (DETWC) arc historically important UBV 
photometric surveys of the Tr 16 and Tr 14 regions. Figure [3] shows the color-magnitude diagram for the deeper 
DETWC survey. It shows that the CCCP X-ray survey detects a roughly constant fraction of the 18 < V < 12 
DETWC stars at all magnitudes and (moderate) absorptions, from K through O stars. It is likely that most other 
CCCP Carina members have similar visual photometric properties. 

Two deep optical UBVRI photometric surveys have been made of a portion of the widely-dispersed Col 228 
cluster around 10'*43'"-60°00' with no strong concentration of OB stars. [Carraro fc Patat| ( |200l] ) (CP) located 1112 
stars in a ~ 20 arcmin^ field, and Delgado et al. (2007) (DAY) located a larger sample in a ~ 100 arcmin^ of which 
152 are identified as likely pre-main sequence cluster members. Unfortunately, this location lies off-axis between two 
ACIS pointings where the X-ray sensitivity is reduced. Only 27 of the CP stars and 17 of the DAY stars are detected 
in X-rays. While small samples, they provide a glimpse at the lower mass function of the stellar population in the 
southwestern part of the Carina complex. 



Table 6 — Continued 



Column Label 


Units 


Description 


J_M(2MASS) 


mag 


J-band magnitude 


H_M(2MASS) 


mag 


H-band magnitude 


K_M(2MASS) 


mag 


Ks-band magnitude 


Ph_Qual(2MASS) 




Photometric quality flag 


CC_Flg(2MASS) 




Contamination and confusion flag 


Id(SOFI) 




row identifier (recno in published catalog) 


Jmag(SOFI) 


mag 


J-band magnitude 


Hmag(SOFI) 


mag 


H-band magnitude 


Ksmag(SOFI) 


mag 


Ks-band magnitude 


Id(NACO) 




row identifier {recno in published catalog) 


Jmag(NACO) 


mag 


J-band magnitude 


Hmag(NACO) 


mag 


H-band magnitude 


Ksmag(NACO) 


mag 


Ks-band magnitude 


Id(Sana) 




row identifier {recno in published catalog) 


Hmag(Sana) 


mag 


H-band magnitude 


Ksmag(Sana) 


mag 


Ks-band magnitude 


Id(SpVela) 




row identifier {DESIG in published catalog) 


mag3_6(SpVela) 


mag 


3.6 /jm magnitude 


mag4_5 (SpVela) 


mag 


4.5 /im magnitude 


mag5_8(SpVela) 


mag 


5.8 /im magnitude 


mag8_0(SpVela) 


mag 


8.0 lira magnitude 


Id(SpSmith) 




row identifier {DESIG in published catalog) 


mag3_6(SpSmith) 


mag 


3.6 /im magnitude 


mag4_5(SpSmith) 


mag 


4.5 lira magnitude 


mag5_8(SpSmith) 


mag 


5.8 lira magnitude 


mag8_0(SpSmith) 


mag 


8.0 lira magnitude 


Id(AC) 




row identifier {CXOTrld in published catalog) 
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4.2. Proper motion measurements of Chandra stars 



The dense, all-sky PPMXL ( [Roeser et ar]|2010[ ) and UCAC3 ( [Zacharias et al.||2010[ ) (with its Bright Star 
Supplement) proper motion surveys have counterparts to 1,006 and 1,471 CCCP sources, respectively. The Carina 
proper motion survey of Cudworth et al. ( 1993[ ) (CMD), based on a century of photographic plates, has 141 CCCP 
counterparts. 

In principle, proper motion measurements can help discriminate rapidly-moving foreground Galactic field stars 
from Carina members. For example, almost all of the PPMXL sources that are spectroscopically confirmed OB 
stars (Skiff 2009| — objects very likely to be Carina members — are found to lie in a small circular region of proper 
motion spac^^that contains a minority (445/1006) of all PPMXL sources identified in the CCCP catalog. However, 
we were unable to identify a boundary in proper motion space that would select a set of CCCP sources that is 
clearly consistent with a foreground population, i.e., an unclustered sample with median X-ray energies indicating 
low interstellar absorption. This is difficult to interpret, and may reflect observational limitations in automatically 
measuring proper motions in nebular regions. As we find no simple interpretation to the proper motion characteristics 



of CCCP sources, we do not consider proper motion to be a useful measurement for classifying CCCP sources (Broos 
eral][20ll| . 



4.3. Chandra counterparts for near-infrared surveys 

Near-infrared (JHK) photometry provides vital information to many CCCP studies. The ubiquitous Two 
Micron All Sky Survey (2MASS) offers very well calibrated positions and photometry across the entire sky, however 
its moderate spatial resolution and sensitivity are not adequate to detect CCCP sources that are crowded and/or 
faint. For the majority of CCCP sources (those in the major clusters) excellent near-IR data from the High Acuity 



A circular region centered at (/^c<,Mi) = (—6.45,-1-2.61) mas yr ^ with a radius of 13 mas yr ^ contains 109 of the 110 PPMXL 
sources that are spectroscopically confirmed OB stars. 




Fig. 3. — BV color-magnitude diagram of V < 18 stars around Tr 14 and Tr 16 reported by [DeGioia-Eastwood et 
al. (2001). All DETWC sources are shown as dots; CCCP detections are shown as pluses. 



14 



Wide-field K-band Imager (HAWK-I) instrumenlj^ ( |Kissler-Patig et"aL 2008 ) was available to us; Preibisch et al. 
(2011 ) discuss these data in depth. Very high resolution, but very narrow-field, observations of the core of the Tr 14 



(2007 


) and by 


Sana et al. 


(2010 



4.4. Chandra counterparts for wide-field mid-infrared surveys 

Photometric observations of mid-infrared (MIR) point sources are invaluable to studies of star forming regions 
in a number of ways. Young stellar objects are readily identified by the MIR emission of their circumstellar disks 
or accreting envelopes, arising from the reprocessing of stellar radiation. Strong constraints on the luminosity of 
OB stars can be obtained from modeling of visual (UBV) and IR spectral energy distributions (SED's), allowing 
the detection of candidate massive stars and the direct measurement of visual extinction ( Povich et al.pOlTal ) . The 
CCCP had access to two data sets from the Spitzer Space Telescope: deep observations of two large fields (Spitzer 
Proposal ID 3420, PI N. Smith) covering the South Pillars and the western wall of the southern superbubble lobe 
( Smith et al.||2010 l, and the shallower Vela-Carina Survey {Spitzer Proposal ID 40791, PI S. Majewski) covering the 
entire field ( [Povich et aL]|2011b[ ). 



5. RELIABILITY OF IDENTIFIED COUNTERPARTS 

Any counterpart matching enterprise is forced to balance two type of failures — missed matches and spurious 
matches — analogous to the balance any source detection enterprise must strike between incompleteness and spurious 



detection. Broos et al. (2007 Appendix) describe a Monte Carlo method for estimating the performance expected 
from our matching procedure. The X-ray catalog is modeled as a mixture of two populations: an "associated 
population" for which true counterparts exist in the counterpart catalog, and an "unassociated population" for 
which no counterparts exist in the counterpart catalogj^ This mixture is represented by the left-hand oval in 
Figure |4] 



X-ray Catalog 



Matching Outcomes 




Fig. 4. — Diagram of all possible ways in which matching an X-ray catalog, comprised of two source populations 
(left oval), with a counterpart catalog can produce associations (upper portion of right oval) and non-associations 
(lower portion of right oval). Section [s] defines the outcome pathways (arrows) and their probabilities (CM, IM, FN, 
FP, TN), and the catalog populations (associated and unassociated) . 

In our Monte Carlo method, the process of matching the associated population to the counterpart catalog is 
simulated many times in order to estimate the frequency of three possible outcomes — correct matches (CM), incorrect 



HAWK-I observations were obtained on the ESO 8-meter Very Large Telescope (VLT) at Paranal Observatory, Chile, under ESO 
programme 60.A-9284(K). 

Some counterparts will be missing from the counterpart catalog, e.g., because they are too bright, too dim, too crowded with other 
sources, or lost in diffuse emission. 
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matches (IM), and false negatives (FN). An IM occurs when the true counterpart is identified as a plausible match, 
but is rejected in favor of an "interloper" — another source in the counterpart catalog that (by chance) produces a 
more plausible match that is mistakenly accepted. A FN is a catalogued counterpart that is in fact associated with 
the X-ray source, but is not identified because the observed X-ray and counterpart positions of the source are too 
distant to meet the match criterion. Separate simulations for the unassociated population estimate the frequency of 
two possible outcomes — true negatives (TN) and false positives (FP) . A FP is a chance superposition of an unrelated 
source in the counterpart catalog and an X-ray source that has no true counterpart. 

The fraction of sources in the actual X-ray catalog expected to belong to the associated population, is 
estimated by equating the number of actual sources for which no match was identified with the sum of our predictions 
for the two pathways that can produce that outcome (the FN fraction multiplied by the size of the associated 
population plus the TN fraction multiplied by the size of the unassociated population); 

[N^at - M) ^ N,at[fA X FN +{1- /a) xTN], (3) 

where N^at is the number of entries in the counterpart catalog and M is the number of X-ray sources that produced 
matches (see Figure |4| . 

We have applied this method to estimate the expected outcomes for the ACIS/Sana, ACIS/2MASS, and 
ACIS/HAWK-I matching procedures for two complementary samples of the ACIS catalog: a set of very signifi- 
cant "primary" sources, defined as having a low probability of being a spurious source (ProbNoSrc_min in Table [l] 
less than 0.003) and a set of less significant "tentative" sources (0.003 < ProbNoSrc_min < 0.01). The very deep 
Sana catalog covers a small field on Tr 14; Figure [5] shows the very good agreement between the CCCP catalog and 
IR observations in this very crowded field. The deep HAWK-I catalog covers about one third of the CCCP field 
(Figure [2|. Although the 2MASS catalog spans the entire CCCP field of view, we estimated matching outcomes over 
only the HAWK-I field of view to facilitate comparison with the HAWK-I results. 

For the primary and tentative samples and for each counterpart catalog. Table [7] reports the total number of 
ACIS sources (column Ncccp) in the field of view analyzed and the expected matching outcomes for those sources: 
correct matches (CM), incorrect matches (IM), and false negatives (FN) for the associated population; true negatives 
(TN) and false positives (FP) for the unassociated population. 

As expected, the fraction of ACIS sources (column CM) expected to produce matches that identify true coun- 
terparts varies with the sensitivity of the counterpart catalog: 81% (primary) and 38% (tentative) for Sana, 74% 
and 23% for HAWK-I, 50% and 12% for 2MASS. Among the primary ACIS sources, the fraction of asserted matches 
expected to be spurious is low: 0.10 for Sana ([6 + 3]/[81 + 6 + 3]), 0.18 for HAWK-I ([9 + 7]/[74 + 9 + 7]), and 0.04 
for 2MASS ([0 -I- 2]/[50 + + 2]). 

Table 7. Estimated outcomes of ACIS/Sana, ACIS/HAWK-I, and ACIS/2MASS matching. 



CCCP Sample Associated Unassociated 



Reliability 


Ncccp 


CM 


IM 


FN 


TN 


FP 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 








Sana 








primary 


380 


81% 


6% 


1% 


8% 


3% 


tentative 


22 


38% 


3% 


0% 


36% 


22% 








HAWK 


-I 






primary 


6749 


74% 


9% 


1% 


9% 


7% 


tentative 


663 


23% 


5% 


0% 


29% 


43% 






2MASS 


in HAWK-I field 




primary 


6749 


50% 


0% 


1% 


48% 


2% 


tentative 


663 


12% 


0% 


0% 


83% 


4% 



Note. — Matching outcomes (Section |5]l consist of correct 
matches (CM), incorrect matches (IM), and false negatives (FN) 
for the associated population; true negatives (TN) and false pos- 
itives (FP) for the unassociated population. 
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Fig. 5. — Matches identified between the CCCP (red diamonds) and a high-resolution observation of the Tr 14 core 
(Sana et al. 2010) (red pluses), with unmatched CCCP (blue diamonds) and Sana (green pluses) sources. CCCP 
sources likely to be outside the Sana field of view (green circle) are shown as black diamonds. The underlying image 
is the CCCP data, scaled so that pixels with a single X-ray event are light gray. The electronic version of the figure 
can reveal, when zoomed, detail not visible in most printed versions. [[Note to referee and journal staff: we would 

like to examine proofs of this figure before deciding whether a second panel, showing a close-up of the most crowded 
region, should be added.]] 



Among the tentative ACTS sources, spurious matches are expected to be more common. Spurious match fractions 
are 0.40 for Sana ([3-h22]/[38-|-3-f 22]), 0.68 for HAWK-I ([5+43]/[23-(-5-|-43]), and 0.25 for 2MASS ([0-^4]/ [12-^0-^4]). 
At least three effects can explain this result. First, the position errors assigned to the tentative ACIS sources (median 
0.44") are larger than for the primary sources (median 0.33"), directly increasing the "footprints" around ACIS 
sources in which a spurious counterpart can be identified. Second, by definition the tentative sample is expected to 
harbor a higher fraction of spurious ACIS detections, for which there is no true counterpart to be found. Third, since 
tentative sources tend to be less X-ray luminous than primary sources, if there is any correlation between X-ray and 
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visual/IR luminosity then tentative X-ray sources are likely to have weak visual/IR emission that is not detectable 
in the counterpart catalog. 

The notion that the X-ray catalog can be viewed as a mixture of associated and unassociated populations — 
on which the matching outcome estimation method above is based — is visually demonstrated in Figure |6j where 
ACIS/HAWK-I separations (in arcseconds) of matching sources are shown for the two simulated populations and for 
the actual ACIS catalog. The upper panel shows that the tentative ACIS catalog (solid) has match separations that 
appear to be a mixture of the associated population (dashed) and the unassociated population (dotted) . The middle 
panel shows the same quantities for the primary ACIS catalog; the fact that the actual ACIS/HAWK-I separations 
(solid) are smaller than those from the associated simulation (dashed) suggests that the adopted ACIS and/or 
HAWK-I position errors are over-estimated. Perhaps a more informative measure of source separation is the smallest 



significance threshold, a parameter to the matching algorithm (see BIO Section 8), required to assert the match. 
The bottom panel presents the distribution of this quantity (expressed as a multiple of a for a Gaussian distribution) 
for the same catalogs and simulations as the middle panel. Separations among matches for the associated population 
(dashed) most commonly fall at ^Ict, falling off at higher separations. In contrast, the occurrence of separations in 
the unassociated population (dotted) increases with separation, reflecting the quadratic growth in area of annular 
footprints with increasing distance from the ACIS source. As in the middle panel, the actual ACIS/HAWK-I matches 
(solid) are generally smaller than those from the associated simulation (dashed), suggesting that the adopted ACIS 
and/or HAWK-I position errors are over-estimated. 



6. QUALITY OF THE X-RAY CATALOG 

6.1. Sensitivity to Point Sources 

We find that very precise language is required to avoid miscommunication when discussing the "sensitivity" of 
a source catalog. Here, we restrict our use of that word to mean the general concept that, for any given observation, 
the probability of detecting a source depends on how "bright" it is; we do not use the word "sensitivity" to refer 
to any specific threshold. We use the term "completeness" to refer to the probability that members of a precisely 
defined parent population of point-like objects will be present in a precisely defined sample of detected sources. In 
contexts where the parent population under discussion is defined by some observational or astrophysical quantity, 
we use the term "completeness limit" to mean a lower-limit on that quantity for which nearly 100% of the parent 
population will be present in a specified sample. Thus, a "completeness limit" is properly stated in terms of a specific 
measured or inferred source property and in terms of a specific sample of detected sources. 

Most studies that involve X-ray-selected sources will have to consider to what extent the analysis is affected by 
the finite sensitivity of the X-ray observations. Sensitivity is perhaps best studied via Monte Carlo simulations of 
synthetic data sets that contain artificial sources that are subjected to the same detection procedure used on the 
observed data. Since most studies need to characterize sensitivity in terms of an astrophysical, not an observational, 
quantity such as intrinsic (corrected for absorption) X-ray luminosity or stellar mass, the first simulation task is 
astrophysical — modeling the X-ray photons that Chandra should see for various astrophysical objects. This task is 
often quite difficult, requiring information that is not well known (e.g., the 3-D distribution of absorbing material and 
the 3-D distribution of sources) and/or requiring astrophysical models that are not well known (e.g., the relationship 
between stellar mass and X-ray luminosity). 

After astrophysical models have predicted apparent (diminished by absorption) X-ray flux for synthetic stars, the 
second simulation task is instrumental — modeling the event data (including background) that would be produced, 
the data reduction process, and the source detection process. If the layout of a Chandra mosaic (exposure times 
and the pattern of observation overlap) is simple and the source detection procedure is simple, then apparent flux 



sensitivity can be studied without detailed Monte Carlo simulations; see for example Georgakakis et al. (20081 and 
references therein. However, for the complex data set and the complex detection procedure employed in the CCCP, 
adequate resources are simply not available to perform the complex simulations that would be required to map 
detection probabilities across the field. 

Instead, we aspire to present here only an empirical description of detection completeness limits with respect 
to apparent photometric quantities, based on the astrophysical assumption that those quantities follow a powerlaw 
distribution that "rolls off" at the low-flux end due to decreasing detection probability. Our primary goal is to 
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Fig. 6. — Observed separations (solid histograms) between tentative (upper panel) and primary (middle and lower 
panels) ACIS sources and HAWK-I counterparts, compared with simulated separations for an associated population 
(dashed) and an unassociated population (dotted). Separations are characterized via angular distance (upper and 
middle panels) and via match significance (lower panel). 



visualize and quantify the large influence that off-axis angle (Theta in Table [IJ has on sensitivity in the CCCP 
catalog (and indeed in many other Chandra catalogs). This relationship arises from two instrumental effects. First, 
the extraction area, and thus any spatially- uniform background component, of a point source with a reasonable 
aperture (e.g., 90% PSF fraction) is ~100 times larger far off-axis (6* = 10') than on-axis. Second, the Chandra 
effective area far off-axis is ^20% lower than on-axis due to mirror vignettingj^ 

The upper panel of Figure [T] shows histograms of the observational quantity NetCounts_t from Table [l] (net 
extracted counts in the total band, 0.5-8 keV) over six disjoint ranges of off-axis angle. The source sample represented 
in Fig.[7]is mapped in FigurejSj we excluded 6877 sources that may exhibit intrinsic spatial variation in the luminosity 



^See Figure 4.5 in the Chandra Proposers' Observatory Guide (http://asc.harvard.edu/proposer/POG/pog_pdf.htmll 
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function (those lying in clusters identified by Feigelson et alT| ( 201 1 )) or that have a poorly-defined off-axis angle (more 
than 20% variation in Theta among the covering observations). 
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Fig. 7. — Histogram of net extracted counts (upper panel) and apparent photon flux (lower panel) for a sample of 
CCCP sources defined in Section [Oj stratified by off-axis angle. 



We choose to interpret the empirical distribution of NetCounts_t in each off-axis angle slice as arising from 
a powerlaw-distributed parent population of astrophysical sources (young stars in Carina, field stars, extragalactic 
objects) that has been incompletely detected at low fluxes. We performed a maximum likelihood fit of the unbinned 
NetCounts_t data in each slice to a powerlaw model with a low-side cutofp^ and we choose to interpret that cutoff 
value as the onset of incompleteness. The best-fit powerlaw exponent and the inferred completeness limit are reported 
as a and in the middle section of Table |8] 

The lower panel of Figure [7] and right section of Table |8] present the same analysis carried out on a simple 
calibrated photometric quantity, an estimate of apparent (not corrected for absorption) photon fiux in the total 
energy band (0.5-8 keV), ^t, photon defined by Equation [ij The source sample depicted is the NetCounts_t sample. 



Maschberger & Kroupa 



I 2009 I discuss the well-known unbiased maximum likelihood estimator for the exponent of an infinite powerlaw 
distribution and they propose a similar unbiased estimator for a truncated powerlaw. We have implemented these estimators using IDL 
in a tool called mLpowerlaw, which is part of our TARA software package ij http: //www, astro .psu. edu/xray/acls/acls_aiialysis .html^ . 



Fig. 8. — CCCP sample used in Figure [7] and to estimate completeness limits in Table |8) these are sources with 
well-defined off-axis angle that are not obviously associated with clusters. 



less 748 sources that have off-nominal exposure time. The best-fit powerlaw exponent and the inferred photon flux 
completeness limit are reported as a and -Ft, photon, lim- 

Those photon flux limits have been converted to intrinsic (corrected for absorption) luminosity limits in the 
total energy band (0.5-8 keV), shown as column Lt,c,iim, by making the following astrophysical assumptions. The 
emitting gas is modeled as the two-temperature thermal plasma that Getman et al. (20101 adopt as typical for low- 
mass pre-main sequence stars with moderate intrinsic luminosity, specifically the model shown in the L^.c — lO^*^" 
entry of their Table 1. The absorbing column is assumed to be Nh = 0.64 x 10^^ cm~^ corresponding to Ay = 4 
mag, which Preibisch et al. ( 2011p find to be typical for Carina stars. The distance to Carina is assumed to be 2300 



Table 8. Completeness Limits for the CCCP Catalog 
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Note. — 

Col. (1): OfT-axis angle range defining the sample 
Col. (2): Median off-axis angle within the sample 
Cols. (3) and (6): Sample size 

Cols. (4) and (7): Estimated exponent of truncated powerlaw model 

Cols. (5), (8), and (9): Estimated truncation threshold (completeness limit) 
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pc. 

We wish to emphasize that the completeness Hmits in Table |8] are a characterization of source samples defined 
only by the CCCP detection process and by off-axis angle. Once a scientific study has imposed additional selection 
criteria — e.g., source classification, the availability of counterpart information, or the availability of individual X-ray 
luminosity estimates — the completeness limit for the resulting sample can be higher. 

The strong variation in detection completeness with off-axis angle (^-^0.8 dex for photon flux) inferred from Tablejs] 
is responsible for the obvious "egg-crate effect" — a greater density of sources on-axis than off-axis — seen when the 



full catalog is plotted on the sky, as in Figure 4 in Townsley et al. (2011a I. A sensible method for suppressing this 



effect from any particular region on the sky (e.g., for a particular star cluster) is to trim the catalog at the photon 
flux limit corresponding to the largest off-axis angle included in that region. Such an approach is recommended for 



spatial analyses involving counting sources e.g., Feigelson et al. (2011). In that study, the "complete" sample of 3,220 



Carina members is defined by log i^t, photon, lim > —5.9 photon cm ^ s "'^ based on the i^t, photon, lim completeness limits 
shown in Table [51 

However, we must emphasize that a catalog so trimmed is not expected to be complete in other astrophysical 
quantities, such as apparent X-ray luminosity (which depends on the temperature of the source) or intrinsic X- 
ray luminosity (which additionally depends on absorption to the source). In fact, applying a cut in NetCounts_t or 
PhotonFlux_t is harmful in such analyses because some of the discarded sources may in fact lie above the completeness 
limit for the quantity of interest. 

Spurious sources have the potential to mask the incompleteness of legitimate detections at low fluxes, thereby 
biasing estimates of completeness limits. For example, in analyses like those in Section |6.1[ spurious sources could 
"fill in" the deficit of detected legitimate sources at low fiuxes in a way that obscures the departure of the detected 
distribution from a powerlaw model. We can roughly assess the degree to which the inferred completeness limits in 
Table |8] may suffer from this issue by estimating an upper limit on the fraction of sources brighter than those limits 
that have a reasonable risk of being spurious. For example, within the innermost off-axis angle slice (6 = 0-3.8') in 
Table [sj ~10% of the sources above the stated PhotonFlux_t completeness limit are "tentative" detections (defined 
by 0.003 < ProbNoSrc_min < 0.01 in Table [T]). Thus, that completeness limit may be somewhat optimistic. All 
other complete samples in that table contain <3% tentative sources. 



6.2. Spurious Sources 

Our detection process is intentionally aggressive, and thus some weak spurious sources should be expected. 
The ultimate strategy for quantifying spurious detections is perhaps Monte Carlo simulation of the detection pro- 
cess executed on synthetic data sets that match the actual observations as closely as possible. Such simulations 
require an estimate of the spatially- variable background in the observations, which is problematic because by defini- 
tion one must decide which photons come from point sources in order to identify which photons are "background" . 
Such simulations must also include artificial point sources because presumably one of the mechanisms for produc- 
ing a spurious candidate in our detection procedure (Section [2]) is imperfect image reconstruction in the wings of 
bright point sources. Such simulations are infeasible due to the complexity of the Carina survey and 
the complexity of our detection procedures. Even if they could be performed, the result of such 
simulations — a characterization of the spurious source population in the full CCCP catalog — would 
be of little practical use because science analyses invariably involve additional sample selections of 
various kinds (e.g., matching to counterpart catalogs, source classification, selection within a region 
of interest on the sky). Propagating simulated spurious sources through sample selection and into 
science analysis would be a very difficult task. 

Perhaps the best method for gauging whether spurious sources may impact a scientific conclusion is to repeat the 
relevant science analysis on a smaller source sample obtained by applying a stricter criterion for source existence, e.g., 
by applying a smaller threshold on the ProbNoSrc_min (Table [I]) statistic that was used to define the CCCP catalog 
(Section [2|. Often, other criteria that define the source sample under study (e.g., involving counterparts, source 
classification, completeness cuts) have already removed most sources whose existence is uncertain. For example, 
in Section |6.1| we found that "tentative" detections are virtually absent from samples that are constructed to be 
complete in apparent photon flux. 
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Although we do not have an accurate estimate for the number of spurious sources in the CCCP 
catalog, a reasonable upper limit for that quantity is the number of X-ray sources for which no 
true counterpart exists in the deep HAWK-I infrared catalog, defined as the "unassociated X-ray 
population" in Section [5[ Recall that the size of the unassociated population and the size of its 
complement, the "associated population", are conceptually distinct from the number of X-ray sources 
without or with identified matches; the unassociated/associated fractions are estimated via detailed 
Monte Carlo simulations. 

The size of the unassociated population is an upper limit, rather than a fair estimate, on the 
number of spurious X-ray detections because not every legitimate X-ray source is expected to be 
detected by HAWK-I. For example, the ^700 extragalactic X-ray sources in the HAWK-I field of view 

1); 



expected to be detected by the CCCP should be very faint in the near infrared ( Getman et al. 



many will be below the completeness limits of HAWKI-I (J~'21, H~'20, K~19, Preibisch et al 



2011 



2011 



)• 



We have estimated the sizes of the unassociated and associated populations, with respect to 
HAWK-I, for the full CCCP catalog and for several smaller hypothetical X-ray catalogs (within the 
HAWK-I field of view) obtained by adopting more conservative requirements on source significance 
(smaller values of the ProbNoSrc.min statistic described in Section [2]). Figure [9] shows that the size 
of the associated population (dark gray bars), which is assumed here to be a lower limit on the size 
of the legitimate X-ray detections, grows as we accept less and less certain X-ray detections. At the 
ProbNoSrc min threshold we chose for the CCCP catalog (0.01), the associated population appears 
to still be growing. The cost of pushing deeper into the data is, of course, a rise in the number of 
spurious sources, assumed here to be bounded by the unassociated population (light gray bars). 

There is no correct decision regarding the tradeoff' between sensitivity and reliability. We generally 
take an aggressive approach to source detection for a variety of reasons. From a near term perspective, 
for studies of diffuse emission we feel that the added contamination arising from failing to identify 
and mask actual point sources is more damaging than the observing area lost to masking spurious 
point source detections. From a long term perspective, we are painfully aware that multiple decades 
may pass before this target is observed again by an X-ray observatory with an angular resolution 
superior to Chandra. We believe that recording tentative X-ray detections in the literature may 
prove valuable to future investigators, who will perhaps have access to observations in other bands 
that are far superior to what is currently available. 
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Fig. 9. — Estimated size of the associated (dark gray) and unassociated (light gray) components of hypothetical CCCP 
catalogs (within the HAWK-I field of view) defined by various thresholds on source significance (the ProbNoSrc_min 
statistic in Section [21). 
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7. SUMMARY 



A careful analysis of the 38 ACIS-I observations (totaling 1.2 Ms) that comprise the Chandra Carina Complex 
Project (CCCP) identifies and extracts 14,369 X-ray point sources; their positions and basic X-ray properties are 
presented in an electronic table. The catalog appears to be complete across most of the field to a limit of ^20 net 
X-ray counts (corresponding to an absorption-corrected total-band luminosity of ^10'^"^ erg s^^ for a 
typical low- mass pre- main sequence star). Near the centers of the ACIS pointings, where the Chandra PSF 
is optimized, the completeness limit improves to ~6 net counts (corresponding to an absorption-corrected 
total-band luminosity of ^10^^^ erg s^ for a typical low-mass pre-main sequence star). A large fraction 
of the detected sources lies beyond these completeness limits, e.g., more than half produced less than 10 net X-ray 
counts. 



Forty sources have high-quality X-ray spectra (more than 500 net counts); most are known ( Gagne et al.||2011 



Naze et al. 


2011 


Parkin et al. 


2011 



2011a 



massive stars. Five suffer from mild photon pile-up in the detector; a new method for pile-up correction is described 
in Appendix [A} 

Counterparts to the majority (69%) of the X-ray sources are identified in a variety of visual and infrared 
catalogs, including newly-available wide-field and deep surveys in the near-IR (HAWK-I on the VLT) and mid-IR 
(Spitzer) bands. Most (>80%) of the X-ray/counterpart source separations are <1"; the reliability of the counterpart 
identifications is discussed. 

The X-ray catalog and counterpart identifications presented here provide the basis for later studies of the young 
stellar population in the Carina complex. Removal of the events associated with these sources allows study of the 
diffuse emission in this starburst region (|Townsley erar]|2011b[). 



CHANDRA CCD PILE-UP RECONSTRUCTION 



As outlined in Section [3] a few CCCP sources are sufficiently strong X-ray emitters to suffer from photon pile- 
— arrival of multiple X-ray photons with a separation in time and space that is too small to allow each to be 
detected as a separate X-ray event. This phenomenon is ubiquitous among CCD-type detectors in X-ray astronomy 
and prevents standard analysis of source ffux, variability, and spectrum ( Allen et al.||1998" Bautz et al.||2000" ). 



Various treatments of pile-up for the Chandra ACIS instrument have been proposed. Broos et al. ( 1998 ) and 
surely other authors point out that pile-up effects in point source spectra can be reduced by simply discarding events 
in the core of the observatory point spread function, where the incident photon flux is highest, and then applying 
appropriate energy-dependent calibration corrections to account for the discarded collecting area. AE supports this 
approach for example, Getman et al. (20051 perform annular extractions for dozens of piled sources in the Orion 
Nebula Cluster. 

The most commonly used pile-up analysis technique involves including pile-up effects in the observatory model 
employed by "forward fitting" spectral analysis packages such as Sherpa and XSPEC. These tools hypothesize an 
astrophysical model of the X-ray source, predict the instrumental spectrum that should be observed by employing 
observatory calibration information, compare the predicted and observed spectra, and adjust parameters of the 



astrophysical model to improve the prediction. Chartas et al. ( 2000 1 describe the derivation of astrophysical models 



for piled sources, using XSPEC to implement an astrophysical model and using a CCD simulator (Townsley et al. 
2002a) to model the response of ACIS, including pile- up effects. 'Davis (2001) modeled pile-up with a parametrized 



integral equation, based on work by Ballet ( 1999) . This model, which we will hereafter refer to as the "Davis model" , 
has been implemented in all the spectral analysis packages in common use by ACIS observers, and has been the 
sole technique applied to most piled ACIS sources for the past decade. Bayesian statistical methods have also been 
studied ( Yu et al!]|2000a|b Kang et ar]|2003 Yu et al.||2004 ), but are not in common usage among observers. 



li 


http : //cxc . harvard. edu/ciao /why /pileup_intro .html 
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A.l. A New Forward-Fitting Algorithm 



Figure [TO] outlines a new variant to the forward modeling approach. Here a non-physical model of the source 
spectrum with many free parameters feeds an incident photon spectrum to the MARX mirror simulatoip"] which 
produces simulated photons by modeling the Chandra-AClS PSF and the observatory dither patterrj^ A model 
of the ACIS CCD (Townsley et al. 2000 2002a|bp simulates (1) the physical interaction of those photons with 
the CCD, (2) the process of reading out charge from the CCD (including charge transfer inefficiencjj^ effects), (3) 
the process of detecting and gradinj^ X-ray events, (4) the event list cleaning process that is applied to actual 
ACIS data, (5) assignment of event positions on the sky to remove the effects of the observatory's dithered pointing, 
and (6) the extraction of events within the aperture chosen by the observer. Since the simulated photons have 
random arrival times, photon pile-up occurs naturally within the simulation in the form of superposition of electron 
charge clouds within individual CCD frames. The non-physical spectral model is iteratively adjusted until the piled 
simulated spectrum is similar to the observed spectrum, both in shape and in overall count rate. When an observed 
piled-up spectrum is reproduced in a satisfactory manner, the simulation is run one final time with pile-up disabled, 
by allowing only a single photon to arrive in each CCD frame. The reconstruction process thus produces two sets of 
simulated events derived from the same set of incident photons; the first exhibits pile-up due to Poisson arrival of 
the photons and resembles the observed spectrum, while the second is free from pile-up. 
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Fig. 10. — Data flow diagram for reconstructing a pile-up-free ACIS spectrum (bottom right) from a piled observed 
ACIS spectrum (top right). No interpretation is attempted on the non-physical model of an incident photon spectrum 
that reproduces the piled observation. 

In principle, an astrophysical spectral model could be used in this fitting process, and then directly 
interpreted scientifically, without a need to simulate an ACIS spectrum free from pile-up. In such 
a scheme, the simulator would replace the response files that a fitting package like XSPEC uses to 
model the detector. However, our CCD simulator is not calibrated to reproduce the response of 
ACIS perfectly (e.g., time- and spatial variation of the ACIS optical blocking filter transmission, 



http: //space .mit . edu/cxc/marx/ 



During an observation, Chandra executes a slow pattern of motion around the nominal target coordinates known as dither fhttp : 



/ /c xc .harvard. edu / ciao/w hy/dither . ht ml ). 



Townsley et al. 



I 2002a I and 



Townsley et al. 



http : / /cxc ■ harvard ■ edu/ciao/why/ct i ■ html 



(2002b I are available in the Physics database of ADS. 
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time- varying charge transfer inefficiency, time- varying amplifier gains). Thus, we choose to fit a non- 
physical spectral model, then simulate an ACIS spectrum free from pile-up, and then fit that spectrum 
to astrophysical models in the normal way. We assume that small imperfections in the calibration of 
the simulation simply distort the input non-physical spectral model, leaving the pile-up phenomenon 
(superposition of photon charge clouds) in the simulation largely unperturbed. Thus, we expect that 
the simulated pile-up free event list is similar to what ACIS itself would have produced if pile-up 
were not present. 

The observed (piled) spectrum is grouped (e.g., so that each group has a SNR of '--'5) and the non- 
physical spectral model is parameterized by an independent photon flux for the energy range covered 
by each group. Within each iteration of the fitting process, the difference between the simulated and 
observed spectrum in each group is used to revise the corresponding photon flux parameter. This 
process is analogous to Lucy-Richardson image reconstruction; the redistribution of X-ray photon 
energy to observed event energy corresponds to the PSF of a telescope blurring the sky onto an 
observed image. The very flexible, non-physical model of the input photon spectrum we employ 
corresponds to the very flexible, non-physical model of the sky employed by Lucy-Richardson image 
reconstruction, which involves a free parameter for every pixel in the image. In both situations, 
the resolution of the model (size of the spectral group, or image pixel size) is usually smaller than 
the resolution of the instrument (energy resolution of ACIS, or PSF of a telescope). In both cases, 
Poisson noise will add features to the observation (spectrum or image) that are narrower than the 
energy resolution of the instrument, and thus cannot be reproduced by the model. This noise will 
tend to be amplified in the model (appearing as "salt-and-pepper" pixels in reconstructed images), 
as it strives to reproduce the noisy observation. However, the noisy photon flux model we derive is 
not interpreted in any way. It merely feeds a final pile-up free simulation, in which the noise in the 
model cannot produce features in the simulated spectrum that are any sharper or larger than the 
Poisson noise in the observation that sculpted the model. 

The Davis model makes no use of any explicit information about the aperture that was used to extract X-ray 
events from the observation of a point source; indeed, that model can be used with diffuse sources. The original 
description of the model ( Davis|2001 1 suggests a very generous point source aperture, and we are not aware of other 



documents that discuss the use of non-standard apertures. One of the attractive features of the pile-up reconstruction 
method outlined here is significant flexibility when choosing an extraction aperture. For example, this method can 
be applied to a piled source in a crowded region extracted with a reduced aperture to avoid a neighbor, as is done for 
two sources in Table[4] A non-standard aperture can be particularly helpful when a source suffers from severe pile-up, 
when the correlation between observed event rate and incident photon rate weakens or even becomes negativej^ Use 
of an annular aperture in such a case discards the severely piled core of the PSF, restoring a positive correlation 
between observed event rate and incident photon rate, which this method requires in order to derive a model for the 
incident photon spectrum. 



A. 2. Application to Carina Sources 

Detailed characterization of the performance of this spectral reconstruction technique will be performed in future 
studies. Ideally, analyses of piled observations would be judged against pile-up-free reference observations, using more 
appropriate targets that do not suffer the time variability expected for the massive stars exhibiting pile-up in the 
CCCP. In this study, we can roughly judge the reasonableness of the spectral reconstruction method by comparing 
its results to those obtained by the standard pile-up analysis technique, the "Davis model" described above. 

For this comparison we choose the two most piled Carina sources from Table |4] that have the large extraction 
apertures required for the Davis model: WR 25 and QZ Car. We assess the consistency between the methods by 
simultaneously fitting the piled and reconstructed spectra in XSPEC 
and then examining the residuals of the fit. The XSPEC model for the piled spectrum includes a Davis model 
component to account for pile-up 



Arnaud|1996 l to a single astrophysical model. 



See Figure 3 in the The Chandra PiSC Guide to Pileup I http://cxc.harvard.edu/ciao/download/doc/pileup_abc.pdf I. 

More precisely, XSPEC compares each of the two spectra to separate instances of a Davis-|-astrophysical model. The astrophysical 
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Figure 11 shows the resuhing piled (red) and pile-up corrected (black) spectra, best fit models, and fit residuals. 
The standard efi^ects of pile-up are seen in the relationship between the corrected and piled spectra: pile-up has 
reduced the overall count rate, but has hardened the spectrum by creating spurious events at high energy by combining 
two or more low energy photons. The most important aspect of these spectra, however, are their model residuals. 
The magnitude of the residuals for the reconstructed spectra and their similarity to the residuals for the piled spectra 
demonstrate the fidelity of the method introduced here, and provide additional confirmation of the Davis model. 
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Fig. 11. — Observed (red) piled spectra and pile-up corrected (black) spectra of the Carina stars WR 25 (top) and 
QZ Car (bottom). The piled and corrected spectra are fit to a common astrophysical model, using the Davis model to 
account for pile-up in the former. Observed spectra are shown as data points with error bars in the upper sub-panels; 
models are shown as histograms; fit residuals are shown in the lower sub-panels. 
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parameters of the two model instances are tied together. The Davis model component for the reconstructed spectrum is effectively 
disabled by setting the frame time parameter to a very small value, and freezing all model parameters. The independent Davis model 
component for the piled spectrum allows the alpha and psffrac parameters to vary, and sets the frame time parameter to the value that 
appropriately accounts for the "FRACEXPO" effect, as described in the The Chandra ABC Guide to Pileup ( |http : //cxc ■ harvard ■ edu/| 
piao/download/doc/pileup_abc.pdf I. 
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